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Abstract. The recent Chandra observation of the radio source at the center of our Galaxy, Sgr A*, puts new 
constraints on its theoretical models. The spectrum is very soft, and the source is rapidly variable. We consider 
different models to explain the observations. We Hnd that the features of the x-ray spectrum can be marginally 
explained with an advection-dominated accretion flow (ADAF) model while it does not well fit the radio spec¬ 
trum. An ADAF with strong winds (ADIOS) model is not favored if we assume that the wind does not radiate. 
Alternatively, we propose a coupled jet plus accretion disk model to explain the observations for Sgr A*. The 
accretion flow is described as an ADAF fed by Bondi-Hoyle accretion of hot plasma in the Galactic Center region. 
A small fraction of the accretion flow is ejected near the black hole, forming a jet after passing through a shock. 
As a result, the electron temperature increases to ^ 2 x 10^^K, which is about 10 times higher than the highest 
temperature attained in the ADAF. The model is self-consistent since the main jet parameters are determined 
by the underlying accretion disk at the inner edge. The emergent spectrum of Sgr A* is the sum of the emission 
from jet and underlying ADAF. The very strong Comptonization of synchrotron emission from the jet dominates 
the bremsstrahlung from the ADAF, therefore, a very short variability timescale is expected and the predicted 
X-ray slope and the radio spectrum is in very good agreement with the observations. 
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1. Introduction 

The energetic radio source Sgr A* located at the center of 
our Galaxy is now widely believed to be the signature of 
a massive black hole with mass M = 2.6 x IO^Mq (Melia 
& Falcke 2001; Haller et al. 1996; Eckart & Genzel 1996; 
Ghez et al. 1998; Reid et al. 1999; Backer & Sramek 1999). 
Its radio spectrum seems to consist of two components, 
with a break around ^ 50 GHz. The spectral dependence 
is oc for jz < 50 GHz, while above this break there 
is a submm bump which is described by oc up to 
^ 10^ GHz followed by a steep cut-off towards the infrared 
(IR) (Zylka et al. 1992; Serabyn et al. 1997; Falcke et al. 
1998). The upper limits from IR (Menten et al. 1997) and 
ROSAT X-ray observations (Predehl & Triimper 1994) in¬ 
dicate that this source is quite dim. 

On the theoretical side, a number of models have been 
proposed in the past years for Sgr A*. Most models are 
based on accretion onto the central massive black hole. 
Possible sources of accretion material include the stel- 
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lar winds emitted by the nearby massive stars and the 
hot interstellar medium. Since in either case the angular 
momentum of the accretion flow should be small, Melia 
(1992; 1994) proposed a spherical accretion model. In this 
model the accretion flow is assumed to free-fall until a 
Keplerian disk is formed within a small “circularization” 
radius. The main contributors to the radio and X-ray 
spectra are synchrotron radiation and bremsstrahlung, 
respectively, from the roughly free-fall flow beyond the 
small disk. However, spherical accretion is likely to be 
an over-simplihcation, since the accretion flow still pos¬ 
sesses some angular momentum. An advection-dominated 
accretion flow (ADAF) model therefore is more dynami¬ 
cally exact in this sense (Narayan et al. 1995; Manmoto 
et al. 1997; Narayan et al. 1998). The most attractive fea¬ 
ture of the ADAF model is its ability to explain the un¬ 
usual low-luminosity of Sgr A* given the relatively abun¬ 
dant accretion material. This is because most of the vis¬ 
cously dissipated energy is stored in the flow and ad- 
vected beyond the event horizon rather than radiated 
away (Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1994, 
1995; Abramowicz et al. 1995; Ghen et al. 1995; Narayan 
et al. 1997; Ghen et al. 1997). In the application to Sgr 
A*, the radio spectrum is produced by the synchrotron 
process in the innermost region of the disk while the X- 
rays are due to bremsstrahlung radiation of the thermal 
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electrons in a large range of radii 10^ - where 

Rs = 2GM/c^ is the Schwarzschild radius. However, the 
ADAF under-predicts the low-frequency radio emission of 
Sgr A* by over an order of magnitude and additional as¬ 
sumptions must be imposed in order to match the spec¬ 
trum (Mahadevan 1998, Ozel et al. 2000). 

Following the initial paper by Reynolds & McKee 
(1980) (see also Blandford & Konigl 1979), Falcke et al. 
(1993) proposed that it is the jet stemming from the disk 
rather than the disk itself which is responsible for the ra¬ 
dio spectrum of Sgr A*. In this model, the submm bump 
is produced by the acceleration zone of the jet, called noz¬ 
zle, while the low-frequency radio spectrum comes from 
the part of the jet beyond the nozzle (Falcke 1996b; Falcke 
& Biermann 1999). The nozzle is of order 10i?s and forms 
from the disk at a radius of ^ 2Rs. This model gives an 
excellent fit to the radio spectrum of Sgr A*, including the 
low-frequency spectrum below the break and the submm 
bump, but the expected X-ray emission was not calculated 
explicitly. 

The latest observational constraints for Sgr A* come 
from the high spatial resolution (w 1 ) Chandra X-ray 
Observatory (Baganoff et al. 2001a, 2001b). Baganoff et al. 
observed Sgr A* twice and they found that Sgr A* comes 
in two states: quiescent and flares. In the present paper 
we concentrate on the quiescent state, whereas the flare 
state is considered in Markoff et al. (2001b). The main ob¬ 
servational results for the quiescent state are summarized 
as follows]^. 

— The absorption-corrected 2-10 keV luminosity is 
(2.2l°;^)xl033erg s-^. 

— The spectrum is well fitted by an absorbed power-law 
model with photon index F = 2.21 q'7. 

— The inner region of the source is rapidly variable on 
short timescale of ~ Ihr. A rapid drop of flux on a 
timescale of 10 minutes is detected in the flare state. 
On the other hand, the comparison between the two 
observations with an interval of about one year indi¬ 
cates that the steady X-ray flux remains almost con¬ 
stant. 

— Some fraction of the X-ray flux may come from a partly 

II 

extended region with diameter « 1 . 

— There is tentative evidence for a Fe Ka line at 6.7 keV. 

These results provide new and strict constraints to the 
theoretical models for Sgr A*. In both the ADAF and 
spherical accretion models mentioned above, the X-ray 
radiation is produced by bremsstrahlung originating from 
10^-10'^i?s- Hence the spectrum is very hard with photon 
index F ^ 1.4 and the predicted variability timescale is 
thousands of hours, much longer than the observed ^ 1 
hour. 

^ The luminosity and especially the photon index are taken 
from Baganoff et al. 2001b, which are slightly different from 
those in Baganoff et al. 2001a where the spectral models used 
did not account for dust scattering; see Baganoff et al. 2001b 
for details. 


Therefore it is necessary to reexamine the theoretical 
models for Sgr A*. Melia et al. (2001) proposed that the 
electrons in the small Keplerian disk can attain a very 
high temperature through some magnetic processes, and 
the resulting synchrotron and self-Compton emission are 
responsible for the radio and X-ray spectrum. However, 
the formation of the small disk may not be a necessary 
result of such low angular momentum accretion. An ac¬ 
cretion flow with very low angular momentum can still 
be described by an ADAF, although such accretion may 
belong to the Bondi-like type rather than disk-like type, 
as shown by Yuan (1999) (see also Abramowicz & Zurek 
1981; Abramowicz 1998). Thus the dynamical scenario of 
this model needs to be studied carefully. 

For the jet model, Falcke & Markoff (2000) take into 
account the contribution from synchrotron self-Compton 
emission (SSC) in the nozzle and find that the parameters 
required to interpret the submm bump give a very good 
fit to the Chandra spectrum without changing the basic 
parameters of the jet model. But the remaining impor¬ 
tant problem in the model is why the parameters of the 
jet possess the required values, particularly in reference to 
the inferred underlying accretion disk. Previous ideas of 
a standard optically thick accretion disk in Sgr A* (e.g., 
Falcke & Heinrich 1994) do not seem to work because the 
predicted IR flux from a standard thin disk with a reason¬ 
able accretion rate would be several orders of magnitude 
higher than the observed IR upper limit (Falcke & Melia 
1997). Therefore, it is crucial to consider the jet and accre¬ 
tion flow as a coupled system in Sgr A*, and to consider 
what are their respective roles if both are truly present in 
Sgr A*. Yuan (2000) presented the first effort, by consid¬ 
ering a combination of jet and ADAF models. However, 
the complete Chandra data was not available at that time 
and the detailed coupling mechanism was lacking in Yuan 
(2000) so it is necessary to revisit the model again. 

The development of the theory provides a new chance 
to model Sgr A*. Since the Bernoulli parameter of the 
ADAF is positive, which means the gas can escape to in¬ 
finity with positive energy, Blandford & Begelman (1999) 
propose an advection-dominated inflow-outflow solution 
(ADIOS) in which most of the gas is lost through winds 
rather than accreted past the horizon of the black hole. 
The concept of strong winds from accretion flow was also 
proposed and studied by Xu & Chen (1997) and Das 
& Chakrabarti (1999). The former described pressure- 
driven winds from centrifugally supported boundary lay¬ 
ers and shocks in the inner regions of disks, and the latter 
proposed an advection-dominated flow where the central 
black hole redirects the inward flow at low latitudes into 
an outflow at high latitudes. We are not explicitly making 
use of the latter two models. The most appealing point of 
the ADIOS model as applied to Sgr A* is that the pre¬ 
dicted X-ray spectrum is possibly much softer than that 
of the ADAF (Quataert & Narayan 1999), and therefore 
could possibly give a better fit to the Chandra data. This is 
because the density profile of the accretion flow becomes 
flatter due to the wind, while X-ray emission at higher 
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frequencies is produced in the inner region of the accre¬ 
tion flow. If we assume that the mass accretion rate in the 
ADIOS is described by a power-law of radius, M oc RP, 
the predicted photon index in Chandra band is approxi¬ 
mately r « 3/2-1- 2p. Thus it is necessary to investigate 
this model for the possibility of interpreting the Chandra 
results. 

In this paper we explore several of the above- 
mentioned models for Sgr A*. By probing a larger pa¬ 
rameter space than before, we find that ADAFs can give 
a marginal interpretation to the new Chandra results, al¬ 
though the fit is not very good in some points (Sect. 2), 
while the ADIOS model can’t (Sect. 3). In Sect. 4 we pro¬ 
pose that the combination of an ADAF and a jet could 
provide an excellent fit to the observations to Sgr A*, and 
present our model results. The last section is a summary 
and discussion. 


2. ADAFs 

We first model Sgr A* with the advection-dominated ac¬ 
cretion model. The modeling technique is described in de¬ 
tail in Yuan et al. (2000; see also Nakamura et al. 1997). 
We use the Paczyhski & Wiita (1980) potential to mimic 
the geometry of the central black hole. A randomly ori¬ 
ented magnetic field is assumed to exist in the accretion 
flow and the ratio between the gas pressure and total pres¬ 
sure (gas pressure plus magnetic pressure) is denoted as 
p. As commonly used, we assume that a fraction, 6, of 
viscous dissipation will directly heat electrons. The ra¬ 
diation mechanisms we consider include bremsstrahlung, 
synchrotron radiation and their Comptonization. We re¬ 
quire a physical global solution that satisfies the no-torque 
condition at the horizon of the black hole, a sonic point 
condition, and the outer boundary conditions. The calcu¬ 
lation of the spectra and the structure of the accretion 
flows are made completely self-consistent as the full set of 
coupled radiation hydrodynamical accretion equations are 
solved numerically. 

The parameters are adopted as follows. We take the 
black hole mass as M = 2.5 x lO^M©, and the viscosity 
parameter is fixed as a = 0.1. We assume the magnetic 
field is in equipartition with the gas pressure or weaker, 
i.e., P = 0.5,0.9,0.99, although the sub-equipartition mag¬ 
netic field is more plausible if a = 0.1. We set 5 as 10“^ 
or 10“^ as usual, i.e., we assume that most of the viscous 
dissipation will heat ions. For the mass accretion rate, us¬ 
ing their latest Chandra observational data, Baganoff et 
al. (2001a) estimate M ^ “i x IO^^Mq yr“^, if the stellar 
wind is the accretion material, or M ~ 1 x 10“®Mq yr“^, if 
the hot ISM around Sgr A* serves as the accretion source, 
which we use as our reference numbers. In principle, the 
accretion rate of Sgr A* could be much higher, based on 
the available material from stellar winds. However, explo¬ 
sive events like the hyper-/supernova Sgr A East could 
temporarily reduce the accretion rate onto Sgr A* sub¬ 
stantially (e.g., Coker 2001). 



Fig. 1. Two fits for a standard ADAF model of Sgr A*. 
The radio and IR data are compiled by Melia & Falcke 
(2001). The short solid line in the X-ray error box shows 
the best fit to the Chandra observation by a power-law 
model in Baganoff et al. (2001b). The parameters for the 
solid line are M = 2.8 x 10“®MEdd, ex = 0.1, P = 0.9,5 = 
10“^. For the dashed line, M = 3 x lO^^MEddjQf = 
0.1, P = 0.99,5 = 10“^. The outer boundary conditions 
are Ti w Te = 8 x 10®K (both), Hout = O.lODxepier (solid 
line) and Dout = O.lbHxepier (dashed line). 

The outer boundary conditions should be taken seri¬ 
ously since they may affect the emergent spectrum signif¬ 
icantly (Yuan et al. 2000). Throughout this paper we set 
the outer boundary of the accretion flow at Rout = lO^Rg, 
where Rs is the Schwarzschild radius of the black hole, 
since this is approximately the location where the accre¬ 
tion begins according to the Bondi-Hoyle capture theory. 
Three outer boundary conditions are the temperatures of 
ions and electrons, Ti^e, and the angular velocity of the 
accretion flow, Dout, at Rout- When Rout is very large, 
as in the present case, the available range of T, e within 
which we can get a physical solution is small, therefore 
the effect of T, e can be neglected. But the feasible range 
of flout is large and may have a significant effect on the 
emergent spectrum. For fixed parameters Rout, P and a, 
there exists a critical value of flout above which the accre¬ 
tion is of disk-like type while below it is of Bondi-like type 
(Yuan 1999). The density of the Bondi-like accretion flow 
is much lower than the disk-like type at the same mass 
accretion rate. Unfortunately the exact value of flout is 
uncertain. We only know that it must be low no mat¬ 
ter whether it originates from stellar winds or from the 
hot ISM. For example, the hydrodynamical simulations of 
Coker & Melia (1997) found flout ~ 0.2 flxepier, if it comes 
from stellar winds. We therefore require in our model that 
flout ^ 0.3 flKepler- 

Fig. 1 shows our fits to the spectrum with an ADAF 
model. The radio and IR data are compiled by Melia 
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& Falcke (2001). The short solid line within the error 
box shows the best fit to the Chandra observation by a 
power-law model in Baganoff et al. (2001b). We show 
rather than uFi, because it is more conducive for judging 
the quality of the spectral fit at radio bands. The solid 
line shows our best fit to the spectrum especially to the 
Chandra data. The parameters are a = 0.1,/3 = 0.9, i5 = 
10-2, and M = 1.55 x IO-^Mq yr"! (or 2.8 x IQ-^MEdd, 
here the Eddington accretion rate is defined as Mfidd = 
TEdd/O.lc^ = 5.525 X IO-^Mq yr”^ where LEdd is the 
Eddington luminosity.). The outer boundary condition at 
.Rout = 10®Rs is Tj Te = 8x 10®K (this value of tempera¬ 
ture is consistent with the X-ray observations by Baganoff 
et al. 2001a), flout = O.lOflKepler- 

Erom the figure we find that an ADAF model can 
fit the X-ray spectrum with a reasonable accretion rate 
^ although the predicted spectrum is flatter than the 
best fit of Baganoff et al. (2001b). The predicted X-ray 
spectrum is composed of two components, namely the 
bremsstrahlung from the outer region of ADAF and the 
second-order SSC from the innermost region of ADAF. 
The dotted line in the figure shows the result exclud¬ 
ing the SSC component. Different from our result, in the 
ADAF model of Narayan et al. (1998), the X-ray emis¬ 
sion is dominated by bremsstrahlung alone. One reason 
for the difference is our use of a higher S, and another rea¬ 
son is that we treat the outer boundary conditions more 
carefully. For bremsstrahlung, the emission at a frequency 
1 / is dominated by the largest radius in an ADAF that 
satisfies hv ~ kT{r). Our numerical calculation results 
indicate that the 2 and 10 keV radiation is dominated 
by radii around 4 x 10'^Rs and 7400Rs, respectively. This 
large radial range is consistent with the extended emis¬ 
sion component (rs 1 w 10®Rs) observed by Chandra. 
In addition, the thermal bremsstrahlung can also explain 
the possible Fe Ka emission line at 6.7 keV (Narayan 
& Raymond 1999). The dynamical timescale of the ac¬ 
cretion flow at these large radii, which is responsible for 
the bremsstrahlung variability, is td « {B?/GMY^'^ ^ 1 
year. Baganoff et al. (2001b) made a comparison be¬ 
tween their two observations with an interval of about 
one year and found that the steady state X-ray flux re¬ 
mains almost constant. This result, combined with the 
rapid variability, seems to indicate that there are two com¬ 
ponents to the X-ray emission operating on very differ¬ 
ent spatial scales and having very different time scales for 
variability. Bremsstrahlung may well be the component 
responsible for the constant flux. The SSC component 
mainly comes from regions 3Rs (see Fig. 1 in Manmoto 


2 In Yuan (2000), the Chandra flux is produced with a higher 
accretion rate, 1.5 x IQ-^MEdd- The discrepancy in accretion 
rate is because in Yuan (2000) the accretion is Bondi-like, while 
in the present paper it is disk-like. For the same accretion rate, 
the density in a Bondi-like accretion flow is much lower than in 
a disk-like flow. Since bremsstrahlung emission is proportional 
to the square of the density the X-ray emission can be largely 
different between the two cases. 


et al. 1997). The corresponding variability timescale is 
^ 3Rs/vr « 1000 seconds. So this component would be 
responsible for the observed rapid variability. 

However, as shown by Fig. 1, this model over-predicts 
the submm bump by a factor of ^ 2-3. We then try 
to lower the synchrotron flux from the ADAF to fit the 
submm bump better, as shown by the dashed line in Fig. 
1. The parameters are M = 1.66 x IQ-^M© yr”^,/? = 
0.99, 6 = 10“®. The outer boundary conditions are the 
same as the solid line except with Hout = 0.15HKepier- 
In this case the second order SSC will become too weak 
to contribute to the X-ray flux, therefore, bremsstrahlung 
is almost the sole contributor to the X-ray spectrum. 
Consequently, the predicted spectrum is too flat and the 
^ 1 hour variability is hard to explain. Considering that 
we can only investigate a limited parameter space of the 
ADAF model, and the fact that the solid line only fits 
marginally, we conclude that it is possible to interpret the 
spectrum of Sgr A* from submm bump to X-ray using the 
ADAF model. However, it remains to be seen whether the 
current ADAF model can indeed produce a strong flare 
as found by Baganoff et al. (2001b). In addition, as in all 
previous ADAF models in the literature, the ADAF model 
always under-predicts the low-frequency radio spectrum 
which needs a contribution from another component such 
as a jet. 

3. ADIOS 

We next attempt to model Sgr A* with an ADIOS. The 
modeling approach is exactly the same as with the ADAF, 
except that the accretion rate is assumed to be described 
by M = Mo(R/Rout)^- We solve the full set of coupled 
radiation hydrodynamical accretion equations to obtain 
the spectra and the structures of the accretion flow con¬ 
sistently. Note that this is an improvement compared to 
Quataert & Narayan (1999) where some dynamical quan¬ 
tities such as radial velocity and sound speed obtained in 
corresponding ADAFs (with M = Mq) are used in calcu¬ 
lating the spectra of the ADIOS. Following Quataert & 
Narayan (1999), we assume that the wind does not radi¬ 
ate. 

We first assume that the fraction of viscous heating of 
electrons is <5 = 10“®. We set a = 0.1 but treat M, p and /3 
as free in order to find the best set of parameters to fit the 
submm bump and the X-ray spectrum. The dashed line in 
Fig. 2 shows our best model results. The parameters are 
Mo = 1.66 X 10-®Mo yr-i, p = 0.28, a = 0.1, but /3 = 0.5 
(not 0.9 since otherwise the predicted radio flux is too low 
compared to the observation). The outer boundary condi¬ 
tions are Ti Ri Te = 8 x 10®K, and Dout = 0.295HKepier at 
Rout = 10® Rs. Compared to the ADAF model, both the 
slope of the X-ray spectrum and the submm bump are now 
fitted better. However, there are two serious problems for 
this fit. The first one is that the required mass accretion 
rate is over 5 times higher than the upper limit estimated 
in Baganoff et al. (2001a) mentioned above. The second 
problem is that the X-ray spectra are produced by thermal 
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Fig. 2. Three ADIOS spectral models for Sgr A*. The 
short-dashed line is for 6 = 10“^,p = 0.28, the long- 
dashed line for 6 = = 0.6. The solid line is for 

6 = l,p= 0.4, the dotted line is exactly the same model as 
the solid line, except that Comptonization of synchrotron 
radiation is neglected. See text for other parameters. 

bremsstrahlung emission alone, therefore this model can¬ 
not explain the short timescale variability. In fact, the in¬ 
troduction of a wind makes the variability timescale even 
longer because the decreasing density of accretion flows 
(e.g. Di Matteo et al. 2000) makes things worse. 

The very rapid variability observed by Chandra indi¬ 
cates that the X-ray emission comes from a very small 
spatial region. This points towards SSC occurring in the 
inner region of the disk. In the case of the existence of 
strong winds, the density of the accretion flow in the in¬ 
nermost region is very low. When the flow is tenuous, SSC 
will show some spectral peaks as a result of different scat¬ 
tering orders. To make SSC dominate over bremsstrahlung 
in the X-ray band, the first order of SSC is more promising 
due to the rapid decrease of Compton scattering probabil¬ 
ity with increasing scattering orders. To make the first or¬ 
der SSC component reach the Chandra band, the electron 
temperature in the emission region must be very high. 

An effective way to increase the electron temperature 
in the accretion flow is to increase <5. In the ADAF we 
generally assume that <5 is as as small as <5 = 10“^ or 
10“^, i.e., the viscous dissipation mainly heats the ions. 
However, because of the uncertainty in the microphysics 
of the ADAF, it is possible that for some reasons, such as 
magnetic reconnection, the viscous dissipation may prefer 
heating electrons, i.e., 5 may be much larger (Bisnovatyi- 
Kogan & Lovelace 1997, 2000; Gruzinov 1998; Quataert 
& Gruzinov 1999; Blackman 1999). In this case, the tem¬ 
perature of the electrons will be greatly increased. 

We try to model the spectrum using various values 
for 5. We find that only when <5 « 1, i.e., almost all of 
the viscous dissipation heats only electrons, can we get 


a high enough electron temperature to make the first or¬ 
der SSC dominate the X-ray emission. The solid line in 
Fig. 2 shows such an example. Other parameters in this 
model are Mq = 2.8 x IO^^Mq yr“^, a = 0.1, (3 = 0.9, and 
p = 0.4. The outer boundary conditions are Ti e = lO^K, 
f^out = 0.25HKepier- The temperature of electrons is as 
high as 10^^ K for the accretion flow within 6i?s and the 
highest temperature is 3 x 10^^ K. This model is then very 
similar to the model proposed by Melia et al. (2001) for 
Sgr A* in the sense that a high-temperature inner disk 
forms, with Te > lO^^K. Synchrotron emission in this 
hottest region produces the submm bump, synchrotron 
self-Compton dominates the X-ray band and gives a very 
soft spectrum. The thermal bremsstrahlung radiation only 
contributes a small part as shown by the dotted line, where 
SSC is neglected. In this case a very short X-ray variability 
timescale can be expected. 

Putting aside the reality of such a high 6, the Ht is 
not satisfactory on the following points: First, it under¬ 
predicts the low-frequency radio spectrum. Second, the 
predicted X-ray slope is much steeper than the best fit 
of Baganoff et al. (2001b). The third problem is that this 
model over-predicts the flux above ~ 100 GHz by a factor 
of 4-6. We cannot get a better fit no matter how we ad¬ 
just the parameters. Because of the strong self-absorption 
of synchrotron emission, the radio spectrum is the result 
of a super-position of blackbody radiation from the differ¬ 
ent parts of the ADAF with different temperatures. So, 
comparing this model with an ADAF (or ADIOS with 
small (5), we can understand that the main reason for the 
over-prediction is its too extreme temperature making the 
flux of the blackbody radiation stronger. Thus we conclude 
that, if we do not consider the possible radiation of winds, 
the ADIOS model is not favored for Sgr A*. 

However, the approximation that the wind does not 
radiate may be an over-simplihcation. For example, the 
part of the wind originating from the supersonic region 
of the accretion disk will possibly be shocked when it is 
ejected out of the disk. Thus it would reach very high 
temperatures and its radiation could not be neglected. In 
this sense, the wind within the sonic radius will present 
itself as radiative, outflowing plasma-i.e., like the plasma 
jets typically observed in AGN. The model would then 
possibly become similar to our jet-disk model presented 
below. 

4. Jet-ADAF model for Sgr A* 

The idea of combining a jet and an ADAF was proposed 
by Falcke (1999) and Donea et al. (1999). Yuan (2000) 
first worked this out in detail and calculated the spectrum 
of the jet-ADAF system for Sgr A* and some nearby el¬ 
liptical galaxies. There is only scant direct observational 
evidence for the existence of a jet in Sgr A*, from the near- 
simultaneous VLBA measurements by Lo et al. (1998). 
They found that the intrinsic source structure at 43 GHz is 
elongated along an essentially north-south direction, with 
an axial ratio of less than 0.3. However, it is interesting to 
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note that the nearby spiral galaxy M81 has a very similar 
radio core and similar unusual polarization features as Sgr 
A* (Bietenholz et al. 2000; Brunthaler et al. 2001). In this 
source, a jet was clearly observed after many VLBI obser¬ 
vations, with the length of the jet being only ~ 400AU at 
43 GHz (Bietenholz, Bartel, & Rupen 2000). If we consider 
M81 to be a scaled-up version of Sgr A*, as suggested by 
their similarity, there could well exist a jet in Sgr A* as 
well. Of course, the jet in Sgr A* would be less powerful 
and hence smaller, making it difficult to detect because 
of the strong scattering of radio waves within the Galaxy. 
More generally, jets seem to be symbiotic with accretion 
disks (Falcke & Biermann 1995; Livio 1999) and they are 
found in basically all kinds of accretion powered systems. 
In this sense the model presented here may be quite gen¬ 
eral. 

The picture of our jet-disk model presented here is as 
follows. The accretion disk is described by an ADAF. In 
the innermost region, r < tq, where parameter ro is the 
jet location, a fraction of the accretion flow is ejected 
out of the disk and forms a jet. Since in our model tq 
is very small (rp w 2i?s, within the sonic point of the 
accretion disk), the radial velocity of the accretion flow 
is supersonic at this small radius (the Mach Number is 
~ 2 — 3). Therefore, when the supersonic accretion flow is 
transferred from the disk into the jet, which is normal to 
the disk, the plasma will be shocked before entering into 
the jet. The shocked gas passes through a nozzle where 
it becomes supersonic. Then it is accelerated along the 
jet axis through the gas pressure gradient force (the grav¬ 
itational force is ignored since its effect is rather small 
in the supersonic regime far away from the black hole) 
and expands sideways with its initial sound speed. Given 
the initial physical states of the plasma at the sonic point 
(top of the nozzle), we can solve for all the quantities as a 
function of distance from the nozzle, and after calculating 
the density and the strength of the magnetic field, we can 
calculate the radiation of the jet (Falcke & Markoff 2000). 

If, however, there exists a possibility that a substan¬ 
tial fraction of the accretion flow can be transferred into 
the jet directly without being shocked (e.g., the accre¬ 
tion flow outside of the sonic point also goes into the jet), 
we could also envisage a mixture of a relatively cold (un¬ 
shocked, ~ 10^°K in the innermost region of ADAF) and 
hot (shocked, ~ lO^^K, see below for this value) electrons 
in the jet. If the energy transfer between particles is in¬ 
efficient, this kind of mixture could last for a long dis¬ 
tance along the jet. For an emission model we can ignore 
this possibility, because the implied radiation should be 
less than the dashed line in Fig. 3 and the contribution 
to the overall spectrum can be neglected. On the other 
hand, such a mixture of hot and cold electrons may be 
needed when considering the circular polarization of Sgr 
A* (Beckert & Falcke 2001). This might increase the cou¬ 
pling constant between jet and disk. 

The exact physics of the nozzle are difficult to model 
since we are at present unclear as to the physical mech¬ 
anism of jet formation. In this paper we treat the nozzle 


only phenomenologically when calculating its spectrum. 
We simply assume that it consists of a series of cylin¬ 
ders with the same electron temperature but linearly de¬ 
creasing density (increasing velocity) from bottom to top. 
The velocity of the gas at the base of the nozzle is as¬ 
sumed to be 1/5 of that at the top of the nozzle where 
it reaches sound speed. The emission is not very sensi¬ 
tive to the exact value of the initial nozzle speed. The 
main radiation mechanisms are synchrotron emission and 
its Comptonization. The parameters describing the jet in¬ 
clude radius and height of the nozzle, ro and zq, electron 
temperature Te, electron number density Ue, the strength 
of the magnetic field B at the top of nozzle, and the angle 
between the jet axis and the line of sight 6. 

All above are free parameters in the original jet model 
(but most of them have obvious physical constraints to 
their range of values, see Falcke & Markoff 2000). But 
here in our coupled jet-disk system, more constraints are 
required so that the jet parameters are consistent with the 
underlying accretion disk. Te is calculated self-consistently 
as follows. When some accretion gas passes through the 
shock and enters into the jet, the ordered kinetic energy in 
the pre-shock gas will be converted into thermal energy in 
the shock front. Neglecting the effects of the magnetic field 
on the shock jump condition (we will check the rational¬ 
ity of this approximation later), we calculate the electron 
temperature of the post-shock plasma by the following 
Rankine-Hugoniot relations, namely the conservation of 
flux of mass, momentum, and energy. Written in the con¬ 
ventional notation, they are 


[pv] = 0, 

(1) 

[P + = 0, 

(2) 

[1 9 2F1 

-v^ + -- =0, 

2 5 p 

(3) 


respectively. To obtain Tg immediately after the shock, we 
still need the ratio between the ion and electron temper¬ 
atures, On the one hand, shock heating may favor the 
ions rather than electrons, as isotropization of the bulk 
flow velocities will give to each species a thermal energy 
proportional to its mass. On the other hand, Goulomb col¬ 
lisions, and maybe collisionless processes also, will bring 
about equilibration between ions and electrons tempera¬ 
tures. Determining the value of ^ is a difficult task (see 
Laming 2000 for a recent review). In the context of su¬ 
pernova remnant shocks, Gargill & Papadopoulos (1988) 
predict ^ = 5 from a numerical simulation, while Laming 
et al. (1996) derive ^ = 20 from their fit to observations. 
We set ^ = 10 in our model. 

For a given shock location tq, we first solve the 
radiation-hydrodynamical equations describing the under¬ 
lying ADAF under selected parameters and outer bound¬ 
ary conditions to obtain the pre-shock physical quantities 
at rg. We then substitute them in the above shock rela¬ 
tions to get the post-shock values. Thus we obtain the 
electron temperature Tg in the nozzle. Therefore Tg is not 
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Fig. 3. The jet-disk spectral model for Sgr A*. The dotted 
line is for the ADAF contribution, the dashed line is for 
the jet emission, and the solid line shows their sum. See 
text for details. 

a free parameter in our model, we can change its value 
only through changing the parameters and outer bound¬ 
ary conditions of the underlying accretion disk. We use the 
same “magnetic parameter” (3 as in the ADAF to describe 
the ratio between the gas pressure and total pressure in 
the nozzle to obtain the value of B if temperature and 
density are known. This means that B is no longer a free 
parameter, either. The density, rig, in the jet follows from 
the coupling constant in the jet-disk symbiosis model 
(Falcke et al. 1993). This is dehned as the ratio between 
mass loss in the jet and accretion rate outside rg. For jets, 
qm is typically a few percent and we require gm ^ 1- 

Our best spectral fit is presented by the solid line 
in Fig. 3. The dashed line denotes the emergent spec¬ 
trum from the jet, and the dotted line is from the un¬ 
derlying disk (ADAF). The solid line is their sum. For 
ADAFs, the parameters are M = 8.8 x 10“^ Mq yr“^, 
a = 0.1, /? = 0.95 and S = 10“^. The outer boundary 
conditions are Ti w Tg = 8 x 10® K, flout = 0.27 flKepier 
at 10® Rs- The parameters for the jet are ro = 1-7 Rg, 

O 

zq = 3.5 ro, qm = 0.5% and 9 = 35 , the “calcu¬ 

lated parameters” are B = 23G, Tg = 2.1 x 10^® K, 
and rzg = 2.4 x 10®. The mass loss rate in the jet is 

Mjet = irTgCsUemp = 4.3 X 10“® Mq yr“^, i.e. 0.5% 
of the accretion rate. For (3 = 0.95 and shock location 
To = 1.7i?s, the Alfven Mach number Ma = vjvAi = 
{v!Cs){‘i'n^pl « 10 > 6, here vai = 7?i/(47rpi)^/^ 
is the pre-shock Alfven speed, the magnetic effects are 
weak in the shock transition condition, so our hydrody¬ 
namic approximation to the shock transition condition, 
Eqs. (l)-(3) above, is justihed (Draine & McKee 1993). 

This model hts the spectrum over the whole range of 
frequencies from radio to the X-ray quite well. The submm 
bump is slightly over-predicted, but it is acceptable con¬ 


sidering the variability of the data in this band (Melia & 
Falcke 2001) and the uncertainty of the model. The low- 
frequency radio emission is mainly contributed by the jet 
outside the nozzle. The contribution from the ADAF is 
rather weak and can be neglected. The submm bump is 
the sum of the synchrotron radiation from both the ADAF 
and the nozzle of the jet. We note that the emission from 
the nozzle is much weaker than in the ADIOS with (5 = 1 
presented in the last section (solid line in Fig. 2) although 
the electron temperatures are both ^ 10^®K. Such a dif¬ 
ference is not surprising considering the much smaller spa¬ 
tial scale of the nozzle, rg = 1.75i?s, while in that case, 
there is a larger radial range with high temperature. In 
this sense, an abrupt increase in the temperature profile 
is necessary to model the spectrum. This is naturally sat¬ 
isfied in our jet-disk model by the formation of a jet. If 
instead the nozzle in our model is replaced by a similar 
high-temperature component such as the inner region of a 
disk, since the temperature prohle of the disk is in general 
smooth, the radial range of this high-temperature compo¬ 
nent would be considerable. In this case, we expect that 
the model would greatly over-predict the submm flux, and 
the low-frequency radio spectrum is still hard to explain, 
as indicated by our calculation for the ADIOS model with 
high (5 in Sect. 3. This is the reason why in Melia et al. 
(2001)’s model the authors require an accretion disk as 
small as ~ 5i?s- 

The X-rays are mainly produced in the nozzle by SSC, 
although bremsstrahlung from the ADAF also contributes 
to a small degree. The predicted X-ray spectrum is the 
sum of the very soft SSC from the nozzle and the relatively 
flatter bremsstrahlung spectrum from the ADAF, which 
is in very good agreement with the Chandra data, almost 
identical to the best fit of Baganoff et al. (2001b). The fit 
is also much better than that of the ADIOS with (5 = 1 
in the last section. In both cases, the X-ray emission is 
the sum of bremsstrahlung and SSC, but in the present 
case, bremsstrahlung produces a much flatter spectrum 
than in the case of an ADIOS due to the absence of a 
strong wind. Because of the contribution of SSC from the 
jet, the variability timescale of X-rays should be short, 
t « rg/fjet ~ 10 minutes. This is consistent with the ~ 
1 hour variability timescale determined in the quiescent 
state and is in excellent agreement with the 600 seconds 
variability timescale detected in the flare state. We show 
that it is the variability of the flux from the nozzle that 
causes the huge-amplitude flare (Markoff et al. 2001b). On 
the other hand, since the bremsstrahlung radiation from 
the ADAF also contributes partly to the X-ray spectrum, 
this could explain the possibly detected extended source 
with ~ 10®i?s(~ 1 )) the 6.7 keV Ko; emission line, and 
steady X-ray flux on ~ one year timescale, as we stated 
in Sect. 2. 

We note that the above nozzle parameters, tempera¬ 
ture, spatial size and density, are very close to the “second 
component” in the model of Beckert (fe Duschl (1997), 
which is also responsible for the submm bump of Sgr 




Feng Yuan et al.: A Jet-ADAF model for Sgr A' 


A*. These parameters seem to be the best ones to fit 
the submm bump. It is interesting that the nozzle with 
these parameters will “evolve” naturally into a jet whose 
emission can well reproduce the low-frequency radio spec¬ 
trum of Sgr A*, and the Comptonization of its synchrotron 
emission can produce a very soft X-ray spectrum which 
can fit the Chandra data excellently. In fact, to make the 
up-scattered submm bump extend to the Chandra band, 
an electron temperature as high as lO^^K is needed. The 
peak frequency of this bump is 10^^ Hz. To up-scatter it 
to the X-ray band, v ^ 10^® Hz, the electron Lorentz fac¬ 
tor must satisfy ss 10^®/I0^^ « 10^. This corresponds 
to a temperature of T « ^^nieC^ ~ lO^^K. This value is 
about 10 times higher than the highest temperature that a 
canonical ADAF can reach in its innermost region, but is 
naturally reached when some fraction of accretion matter 
is shocked^. In addition to a high temperature, the spa¬ 
tial size of the dominant emission medium must be small, 
otherwise the model will over-predict the high-frequency 
radio flux as in ADIOS with high <5 (the solid line in Fig. 
2). This is also easily satisfied in the jet model by requiring 
a small rg. In addition to the above parameters, a trun¬ 
cated (no hard tail) electron energy distribution is also 
required in the model, otherwise the synchrotron emis¬ 
sion will extend above the observed IR flux upper limit. 
Beckert & Duschl (1997) simply assume a mono-energetic 
distribution. In our model, a relativistic thermal distribu¬ 
tion, which is highly peaked at 7 « 3.5-^^, is a natural 
result of shock heating (e.g. Drury 1983) since the Mach 
number is not very large in our case, ~ 2 — 3. 

The mass accretion rate of the ADAF in our model, 
8.8 X 1 O“^M 0 yr“^, is only marginally smaller than 
the lower limit of Baganoff et al. (2001a) estimate of 
1 X 1 O“®M 0 yr“^. If we used a higher accretion rate, 
we would obviously slightly over-predict the flux at the 
submm bump band because of the higher flux of the syn¬ 
chrotron emission from the ADAF. 

There are various ways to further evolve the model. 
One is to introduce global winds in the ADAF. The X- 
ray radiation from the disk would be almost unaffected 
but the radio emission from the disk would be greatly de¬ 
creased because of the great decrease in density close to 
the black hole (Quataert & Narayan 1999). But the wind 
cannot be too strong, otherwise the X-ray spectrum would 
be too soft, as we argued in the case of ADIOS model with 
(5 ~ 1. Another modihcation is to assume that the accre¬ 
tion disk is radiatively truncated within rg, the radius of 
the jet formation (Yuan 2000). The physical reason for the 
truncation is that to form a jet, some amount of energy 
is needed. If we assume this energy comes from the un¬ 
derlying disk, the disk will be left cold within rg because 
of the energy extraction (Blandford & Payne 1982). This 
will greatly suppress the synchrotron emission due to the 


^ We note in this context that Falcke (1996b) and Beckert & 
Duschl (1997) demanded that T « lO^^K based on the Sgr A* 
spectrum alone. 


very sensitive dependence of synchrotron radiation on the 
temperature. 

5. Summary and Discussion 

Recent Chandra X-ray observations put new constraints 
on the theoretical models of Sgr A*. The spectrum is very 
soft, the flux is rapidly variable and the source is extended. 
In this paper we consider three different models to explain 
the observational results of Sgr A*. We find that an ADAF 
model can give a marginally satisfactory interpretation 
to the Chandra spectrum and the rapid X-ray variabil¬ 
ity. But our best fit is still not good for the radio spec¬ 
trum in the sense that it over-predicts the high-frequency 
radio by a factor of 2-3 and significantly under-predicts 
the low-frequency radio. We then consider the possibil¬ 
ity of strong winds from ADAFs, i.e., an ADIOS model. 
If the winds are non-radiative and viscous dissipation in 
the accretion flow mainly heats ions, as generally assumed 
in the literature, this model can fit the spectrum ranging 
from submm bump to X-ray quite well. However, it is hard 
to explain the rapid X-ray variability since in this model 
bremsstrahlung is the sole contributor at X-ray band. If 
we assume that most of the viscous dissipation preferen¬ 
tially heats electrons, a rapidly variable X-ray spectrum 
is expected since in this case the X-ray emission is dom¬ 
inated by SSC. But in this case the model over-predicts 
the radio flux above ~ 100 GHz by a factor of 4-6, and the 
predicted X-ray spectrum is much steeper than the best 
fit of the Chandra observations. 

An excellent fit to all the data including low-frequency 
radio can be obtained with a coupled jet-disk model. In 
this model, the accretion disk is described by an ADAF. In 
the innermost region of the ADAF, ~ 2i?s, some fraction 
<Zm (~ 0.5% if any cold jet component is neglected. See our 
discussion in Sect. 4 for the possibility of a cold jet com¬ 
ponent) of the accretion flow is ejected out of the ADAF 
and transferred into the jet. In this process, a shock occurs 
because the accretion flow is radially supersonic before the 
shock. After the shock the temperature of electrons in the 
nozzle (the base of the jet) reaches about 2 x lO^^K. In this 
case, the synchrotron emission in the nozzle largely dom¬ 
inates the submm bump, and its Comptonization dom¬ 
inates the quiescent X-ray spectrum in Sgr A*. The X- 
ray spectrum is soft and the variability timescale is short. 
Out of the nozzle, the jet gas expands freely outward un¬ 
der the force of the gas pressure gradient of gas pressure. 
Furthermore its self-absorbed synchrotron radiation gives 
an good fit to the low-frequency radio spectrum of Sgr A* 
which is hard to explain in ADAF models. The model is 
completely self-consistent. 

The jet in our model produces a slightly inverted radio 
spectrum, as can be understood from the canonical model 
of Blandford & Konigl (1979), with modihcations as in 
Falcke (1996a). In the absence of a shock acceleration re¬ 
gion in the highly-supersonic outer region of the jet, the 
particles retain the highly-peaked relativistic Maxwellian 
energy distribution which is attained by shock heating oc- 
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curring when the radial supersonic accretion flow is trans¬ 
ferred into the vertical direction. On the other hand, the 
electrons in AGN jets typically seem to have a power- 
law high-energy tail after shock acceleration in jets, since 
the Mach number in jets is very high (Drury 1983). In 
that case, a corresponding optically thin power-law spec¬ 
trum at IR/optical frequencies is generally expected, as 
is seen in many AGN and perhaps even X-ray binary jets 
(e.g., Markoff et al. 2001a). In the case of Sgr A*, the 
absence of an optically-thin power-law indicates that, for 
some unknown reason, such high Mach number shocks do 
not occur. If they would occur under certain conditions, 
we should still see an inverted radio spectrum, but we 
would also expect some kind of hard power-law emission 
at higher frequencies (mid-IR to X-rays). 

In addition to the observations we mention in the 
present paper, there are also constraints to the model 
through the frequency-size relationship obtained from 
VLBI observations (Rogers et al. 1994; Krichbaum et al. 
1998; Lo et al. 1998). The jet-disk model can fit this well 
as shown in Falcke & Markoff (2000). 

We therefore conclude that it is possible to present 
a consistent picture of the emission processes associated 
with the central black hole in our Galaxy by combining 
the three basic astrophysical ingredients that have been 
discussed in recent years: Bondi-Hoyle accretion from the 
immediate environment, optically thin accretion through 
an ADAF, and energy extraction and visible emission by 
a plasma jet. Our jet-ADAF model predicts a closely cor¬ 
related variability among sub-millimeter, IR, and X-ray. 
More broad-band observations and monitoring at various 
wavebands (radio, IR, X-rays) will help to judge whether 
it will be possible to establish a standard model invoking 
those elements for Sgr A* in the near future. For exam¬ 
ple, more precise determination of the IR flux will help to 
further discriminate between the jet-ADAF model and the 
pure ADAF model since the former predicts higher IR flux 
than the latter. This will also be crucial for understanding 
the activity in low-power black holes in general. 
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